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^ . We discuss methods to extract neutrino oscillation parameters based on the di- 
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fS| , rectly observable quantities, without reconstruction of neutrino energy. The dis- 
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tributions of muon energies and production angles are compared to Monte Carlo 
predictions made for a set of different neutrino oscillation parameters. The method 



o . 

CT} , is applied to T2K neutrino beam and tested for a set of MC data samples in order 

to evaluate the statistical error. 

<D I. MOTIVATION 

The traditional method of determination of neutrino oscillation parameters requires re- 
construction of the neutrino energy for each measured event. One assumes the target nucleon 
(N . 

(mass M) to be at rest and bound with mean binding energy e b . The neutrino energy is 

. then reconstructed from the knowledge of charged lepton kinematics: its scattering angle G 

Q\ ' and energy Ef (or momentum kf] 

O" 

F fy (M -e b ) + \ (eg - 2Me b + m*) 
(M-e b )-e f + k f cose 

b : 

The interaction is assumed to be quasi-elastic (QE). The methodology presented in the 
final K2K experiment analysis is based on likelihood functions, which are dependent on 
the estimated neutrino energy reconstruction error However at 1 GeV energies the nu- 
cleon's Fermi motion and non-QE interactions introduce significant and difficult to control 
uncertainty [2]. 

The T2K experiment is a long baseline accelerator muon disappearance experiment held 
in Japan. Its main goal is the precise measurement of Am^ and B23 and determination, 
whether the 613 is nonzero. The expected error of the Am| 3 is smaller, than 4% and for 
the 623 is of the order of 1%. The SuperKamiokande detector would collect about 1600 
Up cc events/year if there were no oscillation [jj, Q]. The statistics will be large enough to 



2 



consider different approaches to extract oscillation parameters and in particular to avoid the 
neutrino energy reconstruction step. Actually, what we observe in SK are mostly muons 
and (sometimes) pions and all the information about oscillation signal is contained in their 
distributions. The problem we address is this: which is optimal strategy (e.g. statistical es- 
timator and data binning) to analyse the data? Our considerations and tests were performed 
for the T2K beam and Super- Kamiokande detection power. 



II. THE METHOD PROPOSAL 



The method is based on the Monte Carlo prediction of muon distribution in the detector. 
The events have been produced by the NuWro neutrino generator 5| . The algorithm of 
oscillation parameter extraction is as follows: 

1. Generate a large number of CC neutrino events with a Monte Carlo generator for the 
given neutrino beam (a typical size of the sample is 1000000 events). 

2. Create reference oscillation samples for a set of different parameters (Ara| 3 ,sin 2 (2@ 2 3)) 
using the muon neutrino survival probability Piy^ — > u^). 

3. Impose detector conditions like muon above Cherenkov threshold, no observed pions. 
Approximate condition for the charged particle being detected is the kinetic energy 
above the Cherenkov threshold. Neutral pions are assumed to be always visible. 

4. Produce independent samples of events of the size corresponding to about 6 years 
of data gathering. The sampling has been performed by a product of two probabili- 
ties: P{ffj, — > Vy) (oscillation signal) and 0.01 (normalisation to the right amount of 
events). This is a way to approximately incorporate the statistical error of the overall 
normalisation. 

Each sample represents a typical set of events as they are expected to be seen in SK. 

5. Compare the experimentally measured muon distribution with the reference samples 
in order to extract the oscillation parameters. Two strategies are considered: 

• x 2 estimator is used (the considered bins are expected to contain at least 10 or 
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20 events): 

X N b - 2 ^ n t { ) 

The index % runs through all bins containing data sets for muon exclusive and 
muon+pion(s) events. The symbol iVj stands for the number of measured events 
in each bin whereas is the expected number of events predicted by the Monte 
Carlo simulation. The normalisation constant Nf, — 2 reflects the fact that our 
fit model has two unknown parameters: sin 2 (2@ 2 3) and Am 2 3 . We go through 
templates produced in step (2) and look for a minimum of the x 2 ■ 

Poisson statistical estimator [6] is used (the considered bins are expected to con- 
tain at least 3 events): 

Ni 



2£ 



rn-Ni + Ni In ■ 
rii 



(3) 



III. PERFORMANCE OF THE METHOD 



The space of oscillation parameters is mapped into a two dimensional lattice on the 
(Am.23 ? sm2 (2023)) plane. The lattice intervals were 5 x 10~ 5 [eV 2 ] in Am| 3 and 5 x 10~ 3 in 
sin 2 (2B23). Physical expectation constraints have been used on the parameter ranges (e.g. 
sin 2 (2e 23 ) < 1, Am 2 3 6 [21,29] x 10" 4 [eV 2 ]). 

Tests have been performed for few values of "true" (Am^ ,sin 2 (2G23)). For each one of 
them 200 MC data samples have been created. 

The histograms have been made for several histogram bin sizes and statistical cuts. The 
range of the bins in muon energy has been fixed to 200-1200 [MeV]. For the x 2 statistics 100, 
80, 75 and 50 [MeV] bins have been used. The binning in cos(B At ) has been varied between 
1/4 and 1/20. The smallest bins were at the verge of Super-K resolution 6j . Alternative 
binning in 8^ has been also performed, there were nine bins between and 90 degrees and one 
bin for the backward scattered muons. This binning has been used to extract information in 
the near 1 [Kt] detector in the Super- Kamiokande experiment [?]. Two statistical cuts were 
imposed on the bins: at least 10 and at least 20 events. Bins with lower expected statistics 
had been discarded. Unfortunately, the x 2 statistics occurred to give unsatisfying results. 
Too few samples were identified with the right oscillation parameters. The Poisson method 
has given better results.The tested bins were 100, 80, 75, 50 [MeV] and 1/5, 1/10, 1/15, 1/20 
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[cos(0 M )]. The best result has been found for 50 [MeV] x 1/5 [cos(0 M )] binning. The figure 
IIHI shows the distribution of results on the reference lattice for (Am^ = 24 x lO -4 ^^ 2 ] 
,sin 2 (2e 23 ) = 0.92) and (Am 2 3 = 26 x 10~ 4 [ey 2 ] ,sin 2 (2e 23 ) = 0.96). Each bin on the plot 
gives the number of data samples identified with one point of the lattice. 



Results for Am 2 23 =2.4x1 0" 3 [eV 2 ], sin 2 (2© 23 )=0.92 Results for Am 2 23 =2.6x1 0" 3 [eV 2 ], sin 2 (2© 23 )=0.98 
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FIG. 1: The figure shows how many of the MC data samples have been identified with each point 
of the reference lattice. 



IV. CONCLUSIONS 

The method gives quite good concentration of the results around the expected value. 
Judging from the plots, without inclusion of the systematic errors, results have 1 a areas 
of about ±0.5 x 10 _4 [eV 2 ] in Am^, and ±0.02 in sin 2 (2023). The error in Am| 3 is close to 
the expected uncertainties of T2K[3|, |4], but the uncertainty of sin 2 (20 23 ) is larger. These 
uncertainties depend strongly on the area, in which the true values of oscillation parameters 
are located. For small values of both sin 2 (2@23) and Am| 3 the concentration of the results 
around true value is rather poor. It seems, that this method is better at extracting the 
squared mass difference then the mixing angle. This happens probably due to the error 
of overall normalisation. The position of the oscillation probability maximum depends on 
Am 2 , whereas its depth depends on the sin 2 (20). Thus the total number of recorded events 
should depend more on the sin 2 (20), then Am 2 . 

More effort must be made to find optimal method of binning. Perhaps using non regular 
bin shapes will improve the method's performance. In general a compromise must be found: 
to have many bins in the region sensitive to oscillation signal but also to keep high statistics. 
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We are aware that the systematic errors will add more uncertainty but its evaluation is 
a separate complicated problem. 
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